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We study quantum transport properties of two-dimensional electron gases under high perpen- 
dicular magnetic fields. For this purpose, we reformulate the high-field expansion, usually done 
in the operatorial language of the guiding-center coordinates, in terms of vortex states within the 
framework of real-time Green functions. These vortex states arise naturally from the consideration 
that the Landau levels quantization can follow directly from the existence of a topological winding 
number. The microscopic computation of the current can then be performed within the Keldysh for- 
malism in a systematic way at finite magnetic fields B (i.e. beyond the semi-classical limit B — oo). 
The formalism allows us to define a general vortex current density as long as the gradient expansion 
theory is applicable. As a result, the total current is expressed in terms of edge contributions only. 
We obtain the first and third lowest order contributions to the current due to Landau-levels mix- 
ing processes, and derive in a transparent way the quantization of the Hall conductance. Finally, 
we point out qualitatively the importance of inhomogeneities of the vortex density to capture the 
dissipative longitudinal transport. 

PACS numbers: 73.43.-f,72.15.Rn,73.50.Jt 



I. INTRODUCTION 

As superfluidity and superconductivity, the integer 
quantum Hall effect is a spectacular manifestation of the 
quantum nature at a macroscopic scale. Interestingly, 
these three phenomena share the same property of dis- 
appearance of the dissipation at low temperatures. In 
the two-dimensional electron gases in the regime of the 
quantum Hall effect, the vanishing of the longitudinal 
magnetoresistance R xx for a range in magnetic fields, is 
accompanied by another remarkable feature which is the 
formation of plateaus of the Hall conductance G xy with 
a robust quantization 1 - in integral multiples of e /h. 

There is now a consensus 2 ^ 3 -^ that the existence of G xy 
plateaus is intimately connected with the localization 
mechanism which causes the vanishing of dissipation. In 
a phcnomcnological picture the localized states that do 
not contribute to the longitudinal transport play the role 
of a reservoir for electrons which allows the chemical po- 
tential to sit between Landau levels e„ = (n + l/2)ftu> c 
(here w c is the cyclotron pulsation). On the other hand, 
delocalizcd states are supposed to exist only within small 
energy bands centered around the energies e n . As long as 
the chemical potential lies within the region of localized 
states, the Hall conductance G xy is expected to remain 
constant and is related to the number of filled bands of 
delocalized states below the chemical potential. The ob- 
servation of wide plateaus implies consequently a large 
density for the localized states with a much more signifi- 
cant spread in energy than for the density of delocalizcd 



states. The density of states with a mobility gap is there- 
fore the essential condition for the quantum Hall effect, 2 

Many different theoretical 

approachea 5 ! 6 . 7 . 8 . 9 . 10 ! 11 . 12 ! 13 ! 14 !",^,!? have been pro- 
posed to explain the quantum Hall effect, which can be 
divided in several classes. The most fruitful treatments, 
such as Laughlin's gauge argument^ and the edge 
state picture ^ 16 ' 17 provide an explanation for the 
quantization of the Hall conductance, but rely on the 
assumption that there exist localized states in the bulk. 
Direct diagrammatic calculations have been developed 
as well. Models with point-like impurities are technically 
manageable^ but require the artificial pinning of the 
chemical potential within gaps in the density of states 
to account for the existence of quantized Hall plateaus. 
Further developments to include finite-range disorder 
within the Green functions diagrammatics were pursued, 
but appear more complicated, since it is needed to 
incorporate more and more diagrams with increasing 
strength of the magnetic fields 

It is worth mentioning that the standard procedure 1 ^ 
that consists in averaging over impurity positions is in 
fact questionable in high magnetic fields. Indeed, it is 
normally justified by the physical assumption of ran- 
domness after successive collision events. However, in- 
stead of the chaotic random walk, the electronic motion 
is relatively regular for high magnetic fields and smooth 
disorder potential. At the classical level, this leads to 
the breakdown of the Drude-Lorentz formula^ Quan- 
tum mechanically, this idea has led to the most simple 
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and illuminating approach to explain the existence of Hall 
plateaus, the guiding-center picture^ ' 10 ^ 1 ' 12 ^ 3 ' 14 ^ 5 ' 21 ' 22 
which is a high-field approximation. Within this "semi- 
classical" description, the quantized cyclotron motion en- 
ergy becomes a constant of motion when the disorder po- 
tential is smooth and varies less than the Landau levels 
separation (in order to have negligible Landau levels mix- 
ing). The reason for this is an ineffective energy exchange 
between the two degrees of freedom, the fast cyclotron 
motion, and the slow guiding center motion. Without 
such an exchange, the strong field localization is then as- 
sociated with the drift motion of the guiding center of 
the cyclotron orbit along contours of constant disorder 
potential. The derealization of electronic states is then 
due to the percolation of equipotential lines, that occurs 
at a single critical energy The localization is thus deter- 
mined in terms of the topology of the equipotential lines, 
and the resulting spectral properties vindicate the mo- 
bility gap argument for the appearance of the quantum 
Hall effect in high fields. 

The computation of quantum transport beyond the 
B = oo limit has however not been achieved in a system- 
atic way^ i 11 ' 12 ^ 4 ' 15 Because the matrix elements of the 
current density relate adjacent Landau levels, the mix- 
ing of Landau levels (which is absent at infinite magnetic 
field) needs thus to be taken into account in a rigorous 
way. In this paper, we reformulate the high field ex- 
pansion in the framework of real-time Green functions 
using a particular set of vortex eigenstates. The frame- 
work developed in our paper, although formally equiv- 
alent to the guiding center approach, allows a straight- 
forward and systematic calculation of quantum transport 
at finite magnetic fields B (i.e. beyond the semi-classical 
limit B = oo), and leads to a transparent derivation of 
both the dissipationless longitudinal conductance and the 
quantized Hall resistance. We also show starting from our 
bulk calculation that the conductance can be expressed 
in terms of edge contributions, thus providing a connec- 
tion between the high-field approaches and the edge state 
picture. This result is found up to the third order in pro- 
cesses mixing the Landau-levels, and may be a general 
result valid at all orders as long as the gradient expansion 
theory is applicable. We interpret this as a fundamental 
manifestation of the robustness of the electronic vortex 
states, which gives rise to strong localization properties 
at large but finite magnetic fields and in long-range dis- 
order. 

The outline of the paper is as follows. In Sec. II, we 
consider the simple problem of a single charged particle 
in a magnetic field. We first analyze within a gauge- 
independent framework the origin of the Landau levels 
energy quantization in terms of the wave-particle dual- 
ity. Such a preliminary discussion is needed because this 
origin is usually obscured by the choice of a gauge and 
the high degeneracy of the energy levels. We derive a 
basis of eigenstates in a way which shows that the en- 
ergy quantization follows only from the single- valucdness 
of the wavefunction, i.e. from a topological condition. In 



Sec. Ill, by using this basis of vortex states, we introduce 
within the framework of the real-time Green functions 
the high-field expansion that allows to treat a smooth 
disorder. In Sec. IV we study thermodynamic properties, 
and especially compute the dependence of the chemical 
potential on the magnetic field, that is important to de- 
termine the width of the Hall plateaus for different tem- 
peratures and disorder broadenings. Quantum transport 
properties in high magnetic fields are investigated in Sec. 
V. We discuss and summarize our work in Section VI. 
Some of the technical details such as the derivation and 
the resolution of the kinetic equation have been collected 
in the appendices. 

II. PRELIMINARIES 

In this Section, we consider a single electron of effec- 
tive mass m* and of charge e = — |e| confined in a (xy) 
two-dimensional plane under a perpendicular magnetic 
field B (pointing in the z direction). The influence of 
an impurity potential is postponed to Section III. The 
hamiltonian has the form 



H = (p - ^A) /2m* = ^-ifiV + ^-Aj /2m* (1) 

where A is the vector potential which is determined from 
the equation 

V x A = B = Bi. (2) 

The vector potential is defined by this way up to the gra- 
dient of an arbitrary function A(r) (i.e. a gauge factor). 
Here, for simplicity, we do not consider the spin degree 
of freedom and disregard therefore the Zeeman splitting 
term. The wave functions 'J and the energy eigenvalues 
E arc derived from the stationary Schrodinger's equation 

H ^ = (3) 

Because it is the vector potential A rather than the 
magnetic field B which appears in Hq , there is some free- 
dom in solving Eq. ([3]). In the next two subsections we 
introduce the two most known ways presented in text- 
books to solve it, before analyzing the question of the 
origin of the quantization of the energy spectrum. 

A. Translation-invariant eigenstates 

The traditional way^ 3 - to solve the eigen-problem is 
to introduce the Landau gauge, for which the vector po- 
tential has the form A = xBy. The hamiltonian is then 
also obviously invariant by translation in the y direction 
so that the wave function is sought under the form 

^(x,y) = f P {x) exp(ipy) 
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with p real. The function f p (x) obeys the differential 
equation 



2m* E 




(4) 



typical of an harmonic oscillator. Here we have intro- 
duced the magnetic length 



h 



(5) 



which is related to the quantum of magnetic flux <&o as 
follows 



$ = 2nl%B=^. 

M 



(6) 



Introducing the new coordinate x' = x — Xq with xq 
-pl%, a solution of ^ is sought under the form 



f p (x) = exp 



2/ 2 



h(x>) 



where h(x') is expanded in power series. The condition 
that f p — > for x — > oo imposes that the function h is a 
polynomial of order n and leads to the well-known result 
of the quantization of the energy spectrum 



E n ,p — fu^c 



1 

n+ 2 



(7) 



where cu c = \e\B/m*c is the cyclotron pulsation. Here 
the energy E is independent of the quantum number p. 
This expresses an important degeneracy of the energy 
levels. Finally, the wave functions are thus labeled by 
the quantum numbers n and p and are written as 



*n, P (a;,y) 



exp 



(x - X ) 2 

211 



{ ^_xo ] 



with m a positive or negative integer, condition required 
by the single-valuedness of the wave function. The func- 
tion fm(r) obeys the one-dimensional differential equa- 
tion 




I ± ( r ^l\ + ( 2m * E 
r dr \ dr J \ h 2 

We search the function /,„ under the form 
The function g m is then given by the equation 



= 0. 



(10) 



(11) 



2\m\ + 1 



2m* 
J? 



d 2 g m 
dr 2 



E - (m + |m| + 1) 



dg n 



1%J dr 



0. 



(12) 



Using the change of variable z = r 2 /2l 2 3l the equation is 
rewritten in the form 



JL-i (m + M + i) 



fjn 



0. 



(13) 



The confluent hypergeometric function $(a,7, z) 
obeys the same kind of equation 



d2 ® , ^ d ^ ^ n 
z—r + (7- z) a<P = 0. 

az z az 



(14) 



There exist two linearly independent solutions to this dif- 
ferential equation $(a, 7, z) and z 1 ~ 1 ^(a~ 7+I, 2—7, z). 
The second solution diverges for z = and is therefore 
not physically acceptable. Furthermore, for z — ► 00 



with H n the Hcrmitc polynomial of degree n. Within this 
basis, the moduli |^ / „ iP (x, y)\ have the property of being 
translation invariant in the y direction, and the eigen- 
states ^n tP (x,y) are not square-intcgrable in the plane. 



B. Rotation-invariant eigenstates 

Instead of the Landau gauge, one can use the symmet- 
rical gauge for which the vector potential isA = Bxr/2. 
Introducing the polar coordinates (r, 9) , one can straight- 
forwardly see that the hamiltonian Hq is invariant by any 
rotation around the origin. Then a convenient form for 
the wave function is 



Th) 
T(a) 

with r(z) the gamma function. Using the fact that the 
product rf m (r) has to be square- integrable, we obtain 
that necessarily 



1 



r(«) 







which leads to a — — I where I is a positive integer. Then 
the confluent hypergeometric functions $(— I, \m\ + l,z) 
are proportional to the generalized Laguerre polynomials 
L| m '(z) of degree I. This results in the energy quantiza- 
tion 



*(r,0) = f m (r) exp(im6») 



(9) 



E = hw c [ I + - [\m\ 



i j = hj c (n + X - J (15) 
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with n = I + (|j7i| + m)/2. The corresponding wavefunc- 
tions are labelled by the quantum numbers m and I 



V2h 



1 71 '■ I 



2/ 2 



Here, the moduli \^ m j(r, 1 



(16) 

are rotation- invariant. In 



contrast with the former basis, the present eigenstates are 
square-integrable in the plane and may thus correspond 
to physical states. The solutions describe states for 
which there is an equally probable distribution of orbits 
whose centers lie on a circle^. 



C. Origin of the energy quantization? 

At this stage, let us compare the origin of the energy 
quantization in the two different derivations. In the first 
case, the quantization originates from the vanishing of 
the wave function far from the origin in the x direction. 
According to the second derivation, two integral numbers 
contribute to the energy quantization: the first integer 
(m) results from the condition of single- valucdncss of the 
wave function, while the second integer (I) stems from 
the condition of square-intcgrability. 

From this comparison, the existence of a physical inter- 
pretation for the energy quantization is not clear. How- 
ever, a clue for an unambiguous origin of this quantiza- 
tion seems to be provided by semi-classical arguments. 
According to the classical equations of motion, a single 
electron confined to a plane and in the presence of a per- 
pendicular magnetic field follows a circular orbit of radius 
R c = v/lu c where v is a constant velocity. The classical 
kinetic energy can thus be written as 

E = -mV = ^-m*(R c uj c ) 2 = Huj c ^- (17) 

with $ c = irR 2 B the magnetic flux enclosed by the tra- 
jectory. Making the correspondence between this classi- 
cal expression and the quantum expression [Eq. (J7]) or 
([15])] of the kinetic energy, we get in the semi-classical 
limit n ^> 1 



(18) 



In this limit, the energy quantization can thus be inter- 
preted as a consequence of the quantization of the flux en- 
closed by the cyclotron trajectory. Furthermore, the flux 
enclosed by the trajectory can be seen in semi-classical 
terms as the phase accumulated during a cycle. We thus 
may wonder how the energy quantization is related to the 
phase of the wavefunction in fully quantum-mechanical 
terms. In subsection III El we shall present a third deriva- 
tion, which has the advantage to show explicitly how the 
energy quantization arises only from a topological condi- 
tion in the point of view of the electron's wave-particle 
duality. 



D. Gauge-invariant formulation 

In order to understand the difference between the two 
bases derived previously and the role of the gauge, it is 
useful to work in a gauge-invariant formalism. An im- 
portant consequence of the presence of a magnetic field 
is that the hamiltonian Hq is always complex whatever 
the gauge chosen for the vector potential. This implies 
the existence of a non-trivial phase for the wave function. 
Let us thus seek ^ in the form 



with / and ip two real functions. Substituting this form in 
the equation ([3]) we get the system of coupled equations 



2tt 







-A/ 



2m*E 



K 2 



-f 



(19) 
(20) 



where the functions ip and / have to be determined to- 
gether with the energy E. Note that in a hydrodynamical 
picture^ of quantum mechanics (see also e.g. Ref. |2(1 ). 
f 2 plays the role of the density of a fictuous fluid and 
Eq. Q19[) is nothing else than the equation of continuity 
for the current probability in the stationary case. 

The explicit expression of the phase <p of the wave func- 
tion depends on the choice of a gauge. Note that how- 
ever the system (|19[) - (|20[) is written in a gauge-invariant 
form. Indeed, if we consider the gauge transformation 
A — > A' = A + VA, the phase transforms as 



ip -> (p 



$0 



(21) 



The gauge is thus reabsorbed additively in the expression 
PTj) of the phase factor. Accordingly, the phase factor 
can be typically divided into two parts 

<P = fgd + <Pgi, 

a part tp gd which is gauge dependent, and a part tp g i 
which is gauge independent. Let us denote the combina- 
tion of the phase gradient and the quantity 27rA/$o by 
the gauge-invariant wavevector K, i.e. 



K 



2tt 



Eqs. p^|) - (P0|) are then rewritten as 



-A/ 



/V-K + 2V/-K = 0, 
0. 



K 2 



2m*E 



f = 



(22) 

(23) 
(24) 



The coupling of the continuity equation (|23|) with the Eq. 
(|24[) was bypassed in the former derivations of the energy 
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eigenvalues by invoking the choice of a gauge and sym- 
metry considerations. However, it is clear here that it is 
possible to get the energy quantization without specify- 
ing a gauge but rather by choosing the vectors K. 

Indeed, for the translation-invariant states, we impose 
that the function / depends only on the variable x. If we 
impose in addition V • K = 0, we obtain from Eq. (f^5|) 
that K || y. Then, the vectors K are fully determined 
from Eq. ([22]) which yields V x K = 2irBz/G> in the ab- 
sence of phase singularities. After integration, we obtain 
that 

K=(„ + ^)y, (25) 

where p is a constant of integration which is obviously 
gauge- independent. Substituting this flow for K in Eq. 
(|2"4"|) . we get the gauge- independent Eq. (0| yielding the 
energy spectrum. 

For the rotation-invariant states, we impose that / de- 
pends only on r and V • K = 0. Eq. (|23]) yields straight- 
forwardly that 

K || z x r. (26) 

In the presence of a phase singularity at the origin, we 
have 

V x K = 27rSz/$ + 27rm(5(r)z, (27) 

with m a positive or negative integer ensuring the singlc- 
valucdncss of the wave function 'J. The form (j2~7| to- 
gether with the condition (|26[) determines also an admis- 
sible flow of K given by 

ttB" 

In the symmetrical gauge, the phase associated with this 
flow is given by cp = (p g i = m9, the gauge-dependent part 
ifgd being simply a constant. 

Therefore, contrary to what could be inferred from the 
derivations III Al and III B[ it is not necessary to choose 
a gauge to write down the Schrodinger's equation ex- 
plicitly. Furthermore, the real difference between the 
translation-invariant states and the rotation-invariant 
states is not the gauge because both kind of states can be 
obtained in any gauge. The real difference is in the choice 
of the flow of K, which is intimately related to the sym- 
metry of the probability density f 2 . The set of quantum 
numbers labelling the cigenstatcs arises directly from this 
choice. For instance, the vector K can be either trans- 
lation invariant or rotation invariant. As a result, there 
are many possibilities to choose other flows of K satis- 
fying the system of equations (f2"3"|) - (f2"4"| . It is somehow 
remarkable that whatever this choice, the energy remains 
quantized as (n + 1/2)Huj c . 



(28) 



center R. At the quantum level, the translation invari- 
ance of the hamiltonian Hq reflects precisely the degen- 
eracy of the kinetic energy with respect to the position 
of the center, while its rotation invariance corresponds to 
the cyclotron motion. Apparently, the former eigenstates 
© or (fTo| do not reflect the symmetry of the cyclotron 
motion around an arbitrary point R so that the consid- 
eration of the classical limit is rather tricky. Since they 
do not correspond to the classical picture of the motion, 
it is difficult to appreciate the wave-particle duality. 

Let us thus try to construct a basis of eigenstates which 
has the right properties. For this purpose, we impose that 
the probability density f 2 has the same symmetry as the 
cyclotron motion, i.e. is a function of |r — R| only. As 
previously, we take the freedom of having V • K = 0. 
Then, we get from Eq. (J23J) that 



K || z x (r-R). 
It is straightforward to see that the flow 



K 



Ir-R 



■kB 
$ 



zx(r-R) 



(29) 



(30) 



is admissible from the point of view of the coupling of 
Eqs. (j2"3|) - (l2"4"l) . Here, we have deliberately located the 
phase singularity at the arbitrary position R (rather than 
at the origin as in III Bp . 

Substituting the form (|5T)|) in the Eq. (|2"4"|) we obtain 
a differential equation for / which is similar to Eq. (jTUJ) 
after introducing the new cylindrical coordinates defined 
by p = |r — R| and = arg (r — R). The function 



fm,n( r ) 



(31) 



with m > is a simple solution of this differential equa- 
tion yielding the energy quantization 



1 

m+- 



(32) 



At this stage, it is not needed to entirely solve the dif- 
ferential equation (|10|) as done in III Bl because we have 
already determined by this way enough states. For the 
same reason, we have restricted oursclf to positive values 
of the integer to. 

Then, considering Eq. (|30[) and taking the symmetrical 
gauge A = B x r/2, we deduce that the phase-gradient 
takes the form 



R 



r-R 



! x R 
2/2 ' 



(33) 



We obtain that the normalized wave functions labelled by 
to and R = (A, Y) are fully expressed in the symmetrical 
gauge as 



E. Vortex eigenstates 

In absence of any potential, the electron follows clas- 
sically a periodic cyclotron motion around an arbitrary 



* m ,R(r) 



Cm 


r R 




V2l B 


Cm 


fz-Z 


77 


\V2l B 



n (r-Rr-2.i-(rxR) 

e miarg(r-R) e IT? 



(34) 
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with C m — 1/V27rm!. z = x + iy, and Z = X + zY". 
For practical convenience, we use henceforth the Dirac 
notation 

*™,R(r) = (r|m,R). 

We now argue that the set of eigenstates \m, R) with 
m restricted to positive values only (m > 0) forms a 
semi-orthogonal basis. The overlap between two states is 
calculated in Appendix A, with the result: 

_ (R-R f ) 2 -2tz-(RxR / ) 

(m,R|m',R) =<We . (35) 



According to (|35|) . the states |m, R) are orthogonal with 
respect to the quantum number m, but not with respect 
to the vortex positions R. However, the overlap is almost 
zero when the distance between two vortex positions ex- 
ceeds a few magnetic lengths. Let us prove that the set 
of vortex states (|34|) is complete. We have 



d 2 R^(r|m,R) (m,R|r') 



d 2 H 



y — 



1 fz'*-Z* 



\21r i v v air 

Recognizing the development into series of the exponential function, we perform first the sum over m to get 



(z-Z 
\V2h 



2 + \z' \ 2 +2\Z\ 2 -2Zz* - 2Z* z' 



d 2 R ( r K R ) ( m > R l r ') = e 

m=0 



2 +\z'\ 2 -2z""2 



r 



e b 



2tt1 2 b 



(36) 



(37) 



The remaining integrations over the vortex coordinates 
R yield straightforwardly 



rf 2 R^](r|m,R)(m,R|r') 



27rZ|<5(r-r'). 



(38) 



Using that S(r — r') = (r|r'), inserting from the left |r) 
and from the right (r'|, integrating over r and r', and 
using the completeness relation for the states |r), we get 
from this expression ([38]) the completeness relation 



-^T X! \m,Il){m,K\ 



2ttI 



(39) 



B 



m=0 



We see that the set of states is complete if we associate 
the area 2nlg with the vortex coordinates X and Y. We 
deduce that the degeneracy of the energy levels is equal 
to (2tt1 2 3 )~ 1 per unit area (i.e. for a finite system it cor- 
responds to <&/<£>() where $ is the total magnetic flux in 
the system), in agreement with the standard derivations 
making use of orthogonal sets of eigenstates. 

We have therefore presented a third way of derivation 
of the eigenvalues and eigenstates for the problem of a 
charged particle in an uniform magnetic field. We find 
that the energy eigenvalues (|32|) are degenerate with re- 
spect to a continuous quantum number, namely the po- 
sitions R, which here correspond to the locations of a 
phase singularity. In the limit of infinite magnetic field 
only, the vortex location R corresponds to the classi- 
cal guiding center coordinates. It is worth noting that, 
contrary to previous standard derivations, we did not re- 
sort to any condition of square-integrability to obtain the 



quantization of the energy levels. Here, the quantization 
of the kinetic energy stems only from the quantization 
of the angular momentum (defined with respect to the 
position R). In terms of the wave-particle duality, the 
mechanism at hand is clearly the interference of the elec- 
tronic wavefunction with itself due to the completion of 
a circular orbit. The eigenstates correspond to vortices 
located at the positions R. Their non-orthogonality with 
respect to R reflects the quantum indeterminacy (2-kI\ 
according to the completeness relation) in the positions 
of the vortices. Eq. (|39|) constitutes an essential relation 
that allows to expand any wave function into a superpo- 
sition of the vortex states |m, R). 

Because one of the quantum numbers (m) results only 
from the singlc-valucdncss of the wave function, i.e. 
stems from a topological constraint, it is sensitive to lo- 
cal perturbations only on the scale of the magnetic length 
around the vortex positions R. This property helps for 
example to capture the interplay of a confining poten- 
tial and the Landau levels in a straightforward manner. 
Within the vortex states basis, it is clear that only the 
states whose vortex positions lie near the edges will be 
affected by the confinement. As a consequence, a devia- 
tion of the spectrum from the Landau levels may be only 
significant near the edges. Similar conclusions concern- 
ing the spectrum of finite systems can be reached after 
thorough calculations for particular models as the Teller 
model (confinement by hard walls) or the model of elastic 
confinement (for a discussion, see e.g. Ref. d). On the 
other hand, the bulk states are not affected by the con- 
fining potential, which justifies the common disregard of 
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the influence of edges when computing thcrmodynamical 
quantities of macroscopic systems. One important ad- 
vantage of the present basis is therefore to allow a sepa- 
rate quantum treatment of the edges and of the bulk. De- 
spite the lack of orthogonality, we shall show in the next 
Section that the basis of vortex states is also very con- 
venient for the treatment of any random potential that 
varies smoothly on the scale of the magnetic length in 
high magnetic fields. 

Finally, we note that the set of states (|34[) has been 
derived previousl y 27 ' 28 in the literature within different 
approaches and contexts, namely the coherent states for- 
malism developed by Glauber— (for a review on the the- 
ory of coherent states see e.g. Ref. l3(il) . Indeed, the states 
can be identified as coherent states of a charged 
particle in a magnetic field and they naturally obey the 
algebr a 29 ! 30 expected for these states (the algebra is char- 
acterized by the identities (|55|) and ([55)) ). A peculiarity 
here is that the vortex states (|3~4")l are both coherent states 
(this property comes out through the degeneracy quan- 
tum number R - the vortex location) and eigenstates 
of the hamiltonian (the eigenvalues are associated with 
the quantum number to) ^ The topological nature of the 
energy quantum number to, that is fundamental to cap- 
ture the quantum transport properties in finite magnetic 
fields as will be developed in the present paper, has seem- 
ingly not been recognized in these approaches! 27 ' 28 The 
set of states (f34|) has already been used for the calculation 
of the partition functional and the Wigner function^ 
We shall present a systematic expansion using this set in 
the framework of the real-time Green functions to study 
both thermodynamic and quantum transport properties 
in high magnetic fields. 



III. EFFECT OF A SMOOTH POTENTIAL ON 
LANDAU LEVELS IN HIGH MAGNETIC FIELDS 

A. Matrix elements of the potential 

In this Section, we consider that the electron feels, in 
addition to the perpendicular magnetic field, the presence 
of a potential V{x, y). The corresponding Hamiltonian is 
therefore 



H = H Q + V(x,y), 



(40) 



where Hq is the kinetic part introduced in the former Sec- 
tion. The matrix elements of the potential V expressed 
in the semi-orthogonal basis |to,R) of eigenstates of Hq 
are given by 



V m , m ,(R,R') = y"d 2 r^(r)^ R (r)* m ,. R -(r).(41) 

We need to arrange the expression of these matrix ele- 
ments by making some manipulations. For this purpose, 
we shall consider that 1^(R) can be expanded in a Taylor 



+00 

V(r)=J2-lr-V u YV(u) 



3=0 ■ 



(42) 



Introducing the relative coordinates d = (R' — R)/2 
and the center of mass coordinates c = (R' + R)/2, we 
first make the change of variable r' = (r— c— idxz)/^/2ls 



V m , m '(R, R') = 2C* m C m '(R|R') Jdx'Jdy' [x' - iy'} m [x' + iy'] m ' V (V2l B r' + c + id x zj e - {x ' 2+v ' 2) (43) 



where (R|R') = (to, R|to, R'), and x' (respectively y') lies in the complex plane on the line defined by z — ~id y j \/ 7 21b 
(z = id x j \P11b)- The potential term is then expanded in Taylor series around the complex point c + id x z 



V 



(^i B r' + c + ! dxz)=y^[r'.V u fy(u) . (44) 



3=0 



Now, if the length scale A characterizing the range over which the potential V varies significantly is typically bigger 
than the magnetic length Ib (this occurs for any nonpathological potentials in sufficiently high magnetic fields), the 
series (|44p converges for all finite x' and y' and represents a function which is analytic throughout the complex plane 
with respect to these two complex variables. Using this analycity property, the contours of integration can be deformed 
to the real axes. Then, noting that 

2r' • V u = (x' - iy') (d Vx + id Vy ) + (x' + iy') (d Vx - id Vy ) , (45) 

with u = (j] x ,Vy)i an d using the binomial theorem, we get after integration that the matrix elements of the potential 
can be written as 
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Then, using the fact that c + jdxz = R + d + idxz and applying once again the Taylor expansion to have the 
potential and its higher order derivatives taken at the point R only, we reorganize the different terms appearing in 
the scries (|46| with respect to the order q of the derivatives of the potential V as 



+00 



V m , m ,(R,R') = (R|R')E^1'( R > R ') 

9 =0 



(47) 



The leading contribution to the series is 



(48) 



For the subsequent calculations (Section V), we shall need also the expression for q = 1 which can be written as 



V mjm' ( R R ') = {£^-W(R)-iS 



(R - R) 



W(R) 



/ 777 / 777 

-Wy {d x V(R) + id y V(R)} S mW + l B Jj {d x V(R) - id y V(R)} 5 m . m>+1 . (49) 



r 



We are now in a position to analyze the peculiarities 
of a smooth potential in the high magnetic field regime. 
Due to the exponential cutoff factor (R|R'), the matrix 
elements V m ,m' (R, R') are non negligible only when R' 
is located at a distance from R which is typically smaller 
than V21b, i-e. the modulous |d| is at best of the order 
of Ib- We deduce that 



V mL' ( R R 



V(R) 



(50) 



where A is the length scale characterizing the spatial vari- 
ations of V around the point R. Therefore, when the 
potential V varies smoothly on the scale of the magnetic 
length Ib, A 3> Ib, we have the inequalities 



kL°L'( r r/ ) » CU RR ') » ^l'( R > R ') 



-(1) 



-(2) 



This shows that we have ordered in the series (|4"Tj) the dif- 
ferent contributions to V m ,m' (R, R') by their magnitude 
of relevance (characterized by q). In the next subsec- 
tion, wc exploit this fact to build a systematic gradient 
expansion within the formalism of Green functions. 

Finally, it is worth noting that the vortex states (j3~4"|) 
are labeled by the continuous index R, and therefore the 
set \m, R) is over-complete in the Hilbert space that has 
a countable basis. Nevertheless, according to the com- 
pleteness relation (|39[) it is possible to expand any arbi- 
trary state or operators in the vortex states representa- 
tion. This is precisely what we have done here for the 
potential. Still, it is in principle necessary to be care- 
ful at this level because the vortex states are obviously 
not linearly independent of one another, and thus the 
expansion of V in terms of these vortex states may be 
not unique. However, clearly, the present expansion (|47[) 
which is only valid for a smooth potential in the high 
magnetic field regime (more precisely the condition is re- 
lated to the criterion of convergence of the series (|4"4"|) ) is 



unambiguous: there is a unique correspondence between 
the matrix elements and the operator V . In the frame- 
work of the coherent states theor y 29 ' 30 the unicity of the 
expansion is provided by the constraint that the expan- 
sion coefficcnts depends analytically upon the continuous 
(complex) quantum number labelling the states. Here, 
we have used explicitly the same condition in the course 
of the derivation (between Eq. (|43|) and (f46|) ) of the ma- 
trix elements of the potential. Namely, the criterion of 
convergence of the Taylor expansion of the potential has 
implied the analycity of the function V, what renders the 
expansion unique. 



B. Green functions formalism 

To investigate thermodynamic and transport proper- 
ties of the two-dimensional electron gas in a smooth po- 
tential V and in high magnetic fields, we shall use the 
formalism of real-time Green functions. Retarded, ad- 
vanced and Keldysh Green functions are respectively de- 
fined as (see e.g. Refs. H3 and HI) 

G R {x 1 ,x 2 ) = -i6(h-t 2 ) ({^(x 1 ),^(x 2 )}) (51) 

G A (x 1 ,x 2 ) = i9(t 2 -h) ({*0n),ft(x 2 )}) (52) 

G K {x 1 ,x 2 ) = -i ([tf(a:i), (a*)]) (53) 

where [ , ] means the commutator, { , } the anticommu- 
tator, and 9 is the Heavisidc function. The Green func- 
tions relate the field operator ^(x) of the particle at one 
point xi = (ri,ii) in space-time to the conjugate field 
operator ^^(x^ at another point x 2 = (r 2 , t 2 ). The field 
operators ^(a^) and ^^(x 2 ) are expressed in terms of the 
eigenfunctions ^(r) and eigenvalues 



(54) 
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*t( X3 ) = ^ c t**(r 2 )e^ 



(55) 



where c\ and c„ are respectively the creation and de- 
struction operators. 

In the absence of the potential V, the Green functions 
can be written explicitly. We denote them by Go- Due to 
its relevance for the high magnetic field expansion theory, 
we shall use the basis of eigenstates \v) = \m, R) derived 
in Section II. E, which is complete with the prescription 
(see the closure relation ([59"]) ) 



Ex / gPR 



(56) 



v m— ' 

The different Green functions Go are thus given by 
G§(x u x 2 ) = -^(t)^* 1/ (r 1 )*:(r 2 )e- ie " i , (57) 

a*) = ^(-t)^^(r 1 )*:(r 2 )e- ie " i , (58) 

V 

G${xi,x 2 ) = *^(2n F (e,)-l)4',(r 1 )**(r 2 ) e -^ t 

(59) 

where t = t\ — t% and 

1 



1 + exp(£„/T) 



(60) 



is the Fermi-Dirac distribution function with £ v = e v — /i 
and s v = (m + l/2)hw c . 

Furthermore, instead of working in the representation 
position |r), we shall henceforth work in the representa- 
tion | v) where the Green functions take a simpler form 
due to the semi-orthogonality property of the basis v 



G^{y x ,t x \v 2 ,t 2 ) 
Go(isi,t 1 ;v 2 ,t 2 ) 



-i0(t) (fih)e- ie "i', 
i0(-t) (^|i^>e- fe -x*, 



i (2n F (£„ 



1) (vi\v 2 )e~ 



(61) 
(62) 
i (63) 



The Green functions are diagonal with respect to the 
energy quantum number m but not with respect to the 
degeneracy quantum number R. After Fourier transfor- 
mation with respect to the time difference t, the Green 
functions are written in the energy (uj) representation as 



G^{vi;v 2 ) = 



S mi ,m,2 (R1JR2) 
w - C, ni + iS 

S mi ,m 2 (R1IR2) 

W - C mi - iS 



(64) 
(65) 



Gq{v x \v 2 ) = -2wrtanhQ^ (vi\i/ 2 }5(u - £ mi ). 

(66) 

For many calculations with the basis \v), it is useful to 
note that within each eigensubspace, we have the follow- 
ing identity 



d 2 R 



,R'|TO,R)(m,R|m,R") = (m,R'\m,R" 



(67) 

which can be interpreted as a closure relation that holds 
in the fixed m eigensubspace 



J 2^ |m ' R)(m ' 



Rl = 1. 



(68) 



High magnetic fields gradient expansion: 
retarded and advanced Green functions 



The retarded and advanced Green functions in the 
presence of the perturbation V are obtained from the 
Dyson equation, which takes the form in the energy (u>) 
representation (for more details see Appendix B.l) 



(u-Z mi ±iS)G% A V2 (cj) 



\V\\Vl) 



E 

"3 



(69) 

Note that the quantum-kinetic equation obeyed by the 
Keldysh Green function G K in the representation \v) is 
derived in the Appendix B.l. The function G impor- 
tant for quantum transport properties is determined in 
the Appendix B.2. 

Here V includes both the scalar electrostatic poten- 
tial and the impurity potential. The explicit form of 

V does not matter at this step. We only assume that 

V is a smooth potential in high magnetic fields. As 
mentioned above, the potential term V VlU2 is thus princi- 
pally described by the first terms in the series expansion 
P?)) . The retarded and advanced Green functions are 
also sought as a series expansion of terms of the order of 



+00 

3=0 



(70) 



Inserting the leading contribution to V VlV2 in the Dyson 
equation (f6"9"]) . the zero order Green functions G^ (in the 
presence of the potential) are given by 



(71) 
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with Vu^ u . 3 — V(Ri) S mi _ ril3 . The summation over 1/3 is readily performed using the relation (|68[) so that we get 

ri(o)R,Af, ,\ — ^mi,m 2 (R1IR2) (70 \ 
(j- £ mi - V Ri ) ± id 



We note that in the denominator of the Green functions (|72|) . the potential can be taken cither at the point R x or at 
the point R2 since, due to the presence of the overlap (R1IR2) in the numerator, V(Ri) differs from V^(R2) at best by 
a correcting term of the order of VI b/ A which can be disregarded in the zero order expression of the Green function. 
The absence at leading order of transitions between the quantum numbers m is a clear manifestation of their purely 
topological (and local) nature: they are immune to small perturbations. The leading effect of the perturbation V is 
that the vortex state acquires a local potential energy V^(R) which is responsible for a lifting of the Landau levels 
degeneracy with respect to the quantum number R. 

The next order Green functions G^ in the presence of the potential are given by 

(w - C mi - V(Ri) ± iS) G« f V) = (Ri|R 3 }GWf ». (73) 

Inserting the expression (fT2)) for G < - 0%> with the potential V in the denominator taken at the point R2, and performing 
the summation over R 3 in the right-hand side of ([73)) by considering that 

J2 J 



J ^R 3 (Ri|R3)(R3|R2)-i[R 2 + Ri-^x(R 2 -R 1 )](R 1 |R 2 ), (74) 



we obtain 



r (l)R,A( \_ V r m 1 i ) ,m 2 (Rl,R2) (R1IR2) 

1/11/2 K ' (u-Ui-V(Ri)±i5)(w-U 2 -V(R 2 )±i6y [ ' 



IV. DENSITY OF STATES AND CHEMICAL 
POTENTIAL 

In this Section, we consider an equilibrium situation 
where the potential energy V(R) consists of an impu- 
rity potential only. Our goal is to investigate thermody- 
namical features such as the dependence of the chemi- 
cal potential on the magnetic field under the regime of 
a smooth potential in high magnetc fields. For this pur- 
pose, we first derive the energy dependence of the density 
of states. 



A. Density of states 

The density of states g(e) at the energy e is obtained 
from 



7T ' " ' 



(76) 



where e = lu + [i. Inserting the zero order Green function, 
the total density of states is straightforwardly given at 
leading order of the expansion theory by 

+°° p j2tj 

f( £ )=E /2-3-<y(e-(m+l/2)^ c -V(R)). (77) 

Note that the next correction term to the density of states 
that we do not give here (although it can be derived in 



a systematic way rather easily from our formalism) is 
provided by the second-order Green function G^ . In- 
troducing 



P(V) 



J^6(V-V(R)) 



(78) 



where S is the total area of the system, the density of 
states takes the form 



S 



+00 . 

= TV E / dV P( V ) He-(m+ l/2)tuv c - V) 

77 B m= QJ 



(79) 

The distribution of the amplitudes V depends on the 
model considered. If we assume a gaussian distribution 
to account for the effect of a disordered potential present 
in the bulk 

p{V ) = ^L-e- v2 ^\ (80) 



2tiT 

the density of states corresponds obviously to a sum of 
gaussian peaks 



9(e) 



S 



+00 ^ 

E ~JK, 



2lTl B V27IT 

m— f) v 



(81) 



m=0 

Using the Poisson summation formula 

+00 +00 /*+C 



£/(") = E / f(t)e- 27 ™dt, (82) 

m=0 i=-oo 
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g{£F)/go 




crucial to understand the width of the plateaus of the 
Hall conductance in the integer quantum Hall effect, see 
Section |Vl) . For this purpose we follow the derivation of 
Refs. l38l and l39l by calculating first the grand canonical 
thermodynamic potential f2 



fi 



-2T 



+ 00 



de g(e) In ( 1 + e 



(84) 



Here the factor 2 accounts for the spin. Inserting the 
expression (|83[) for the density of states, and performing 
the integration over the energy by using an integration 
by parts (for details see Ref. I38T ) , we get that at temper- 
atures T <C /i 



FIG. 1: (color online) Density of states at ef as a function 
of the inverse magnetic field for various disorder broadening 
parameters F The solid curve is obtained for F = 0.1 ef, 
the dashed curve for F = 0.03 ef, and the dotted curve for 
r = 0.01 ef. 



(85) 



where the field non-oscillating part of the thermodynam- 
ical potential is 



n 



(86) 



we can express the density of states (per spin) in the form anc j ^he oscillating part is 



5(e) = 5o < 1 + 



2£(-l)< 



i=i 



f2irle\ 



nrtj z — ' 1 v 



(83) 

where .go = Sm* /2irh 2 is the density of states (per spin) 
at zero magnetic field. For illustration, we have plotted 
in Fig. Q]the density of states obtained for three different 
broadening parameters T. 

This result (j83|) for the density of states in high mag- 
netic fields and in the presence of a random potential 
with a large correlation length has been obtained within 
other theoretical approaches, such as the conventional 
perturbation theory— the quasiclassical path integral 
formalis m 36 ' 37 or the high-field expansion^ It should be 
noted that the obtention of a gaussian form for the Lan- 
dau levels broadening in the conventional diagrammatic 
technique is highly nontrivial (and obtained at the price 
of technical difficulties, see Ref. Ii~8h since this technique 
implicitly assumes a Lorentzian broadening of the Lan- 
dau levels. On the other hand, the Gaussian broadening 
is naturally derived from the high-field expansions^ and 
the quasi-classical path integral technique i 36 ' 37 



B. Oscillations of the chemical potential 

One great advantage of formulating the high-field ex- 
pansion in terms of real-time Green functions is the in- 
sight into the quantum mechanical transport properties 
in high magnetic fields (see next Section). To study the 
dependence of the transport coefficients on the magnetic 
field, it is important to determine first the dependence 
of the chemical potential itself on the magnetic field (the 
location of the chemical potential with respect to the 
Landau levels peaks when varying the magnetic field is 



, (87) 

with A; = 2ir 2 lT/hLu c . The number of electrons is ex- 
pressed as 



TV 



on 

<9/i 



T.B 



which yields 



TV 



2.90 > 

TT * 



(-1)' A 



l=i 



I sinh A, 



sm | — | e v I 



+2.90/U- 



(89) 



This equation determines the number of particles TV as 
a function of (B, T) at fixed /i. On the other hand, we 
can consider it in the thermodynamical limit as the equa- 
tion for the chemical potential as the function of (B,T) 
at a fixed number of particles TV. Introducing the Fermi 
energy ep — N/2g , the dependence of the chemical po- 
tential /i on the temperature T and the magnetic field B 
has then to be found from the self-consistent equation 



hijj. 



£ 

1=1 



-iy 



i 



A 



2-kI^l 

huj c J sinh A; 



(90) 

Here the parameter T, which is independent of the 
magnetic field, characterizes the global variations of the 
potential amplitudes in the bulk. For systems with di- 
mensions exceeding the magnetic length I b by several or- 
der of magnitudes the ratio T/huj c can in principle take 
values bigger than unity if the magnetic field is not too 
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(jl/ef 




FIG. 2: (color online) Oscillations of the chemical potential 
as a function of the inverse magnetic field (through the ratio 
EF/huJc) for various disorder broadening parameters V (same 
values as in Fig. [T] we have used the same convention for 
the correspondence between the curves and the parameters). 
The plot is made for the temperature T = £f/200. For large 
broadening, the oscillations are damped. 



high or/and if the system is sufficiently disordered (so 
that the spread of the potential amplitudes exceeds the 
cyclotron energy). In this case, there is a significant den- 
sity between the Landau levels and we see from Eq. (|90|) 
that the chemical potential (i depends weakly on the mag- 
netic field as 



/i ps e F 



huj c 



Ai 



2-KEl 



hojc J sinh Ai 



(91) 



In the opposite regime T <C huj c , the density of states 
consists of very sharp bands peaked on the Landau lev- 
els, and Eq. ([50)) yields that at low temperatures the 
chemical potential is pinned to a value close to e n for all 
magnetic fields except for the integral values of Ef/Hlo c 
where it experiences a sharp jump. 

To confirm these trends at a more quantitative level, 
we have solved numerically the self-consistent equation 
(|9TJ|) . In Fig. we have represented the oscillations of the 
chemical potential as a function of the inverse magnetic 
field for three different values of the disorder parameter 
r (for convenience, we have used the same values for F 
as in the plot of the density of states, Fig. [1]). The char- 
acteristic sawtooth variation (dotted and dashed curves) 
is obtained for a vanishing density of states between the 
Landau levels. On the other hand, the oscillations of the 
chemical potential /i around its zero field value (£ F ) are 
smoothed and damped in the case of a significant overlap 
between the Landau levels. 



V. QUANTUM TRANSPORT PROPERTIES 

A. Derivation of the current 

We consider the problem of electrons confined in a Hall 
bar of finite width L y . The edges are located at the posi- 
tions Y = ±L y /2 and are included in our formalism by a 
rigid confinement brought about by infinitely high walls. 
In experimental devices, a current I x is applied along the 
longitudinal x direction and induces a transverse voltage 
difference Uh (Hall voltage). The hamiltonian consid- 
ered for the bulk of the sample is given by (j4"0)) where 
the potential term V(r) includes both the impurity po- 
tential and the electrostatic potential (we remind that 
the explicit form of the potential is not needed: we just 
assume that V varies smoothly on the scale of the mag- 
netic length). 

We determine the total current I from the bulk Green 
functions. It is expressed as 



1 = 



E 

v.v' 



2^< 



,M (z/|j>) 



(92) 



where L x is the longitudinal length of the sample. The 
matrix elements of the current density operator j are 
given by (see Appendix C) 

_ ie?t(R/|R) ^ y/m8 m , m '+i - VrrJ 8 m ',m+i u^) 



and 



G< = i (G K - {G R - G A )) 



(94) 



We have established in Appendix B that the identity 

G K = tanh(^/2T)(G fi - G A ) (95) 

that relates the Keldysh Green function to the retarded 
and advanced Green functions at equilibrium is also valid 
at any order of the expansion theory in the stationary 
regime. This result corresponds to the realization of a 
nonequilibrium regime with a local (hydrodynamic) equi- 
librium. Therefore 



G< = -n F ((j) (G R - G A ) 



(96) 



It is straightforward to see that the leading order Green 
functions G^ which are diagonal with respect to the 
quantum number m do not contribute to the current since 
this latter involves only terms that couple the quantum 
numbers m' and mil. As emphasized in the Intro- 
duction, it is thus necessary to calculate the first con- 
tribution introducing Landau levels mixing in order to 
determine the quantum transport properties, i.e. it is 
required to work beyond the semi-classical limit B = oo 
(unlike spectral properties that can be already obtained 
in the absence of Landau levels mixing). Then, inserting 
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the first order Green functions G^ , and performing the 
sums over to, R' and m', we get the simple expression 

+00 „ ,2f> 

I= IJ2 — n F (^ m + V(R)) Vy(R)xz. (97) 

Before deriving the quantization of the Hall conduc- 
tance, let us comment on this expression. We see that 
a current exists only in the presence of gradients of the 
potential V. By disregarding the quantization of the cy- 
clotron motion and therefore the sum over Landau lev- 
els, we obtain a formula for the current that is similar to 
the formula that would be obtained classically by solving 
the Boltzmann equation together with the equations of 
motions in the limit of infinite magnetic fields (see Ref. 
ITih . In the limit B — ► 00 the current is described only 
in terms of the drift of the center of the cyclotron or- 
bit (in this limit the positions of the center and of the 
electron are in fact the same since the cyclotron radius 
vanishes). In our quantum-mechanical picture, a current 
is obviously associated with a drift of the vortex location 
only. The contribution of the orbital currents (induced 
by the magnetic field) to the total current is taken into 
account through the sum over the integral quantum num- 
bers m and the quantization of the spectrum into Landau 
levels. Due to this quantization of the phase circulation 
around the vortex that confers in a way a quantum ro- 
bustness to the cyclotron motion, the above current for- 
mula holds at finite magnetic fields in any potential that 
varies smoothly at the scale of the magnetic length. We 
propose now a transparent derivation of the quantization 
of the Hall conductance that differs from that presented 
in the Ref. (although the starting formula (|9"T|) is sim- 
ilar). 

The crucial and remarkable point is that from the for- 
mula (|97p we can formally extract a vortex current den- 
sity j(R) (we remind that here R corresponds to the 
vortex position, and not to the electron position) 

+00 

j(R) = I n F (&» + y ( R )) W ( R ) x ( 98 ) 

m=0 

which can be rewritten as 

m=0 

where 

N F (e) = e + Thxn F (e) (100) 

is an antiderivative of the Fcrmi-Dirac function n F {e). 
We can easily check that the current density (f9"9"| obeys 
the continuity equation V • j = 0. 

The current flowing in the longitudinal direction is 
then given by 

/+00 
dY d Y N F {U + V(R)) . (101) 

m=0 



We can readily integrate this expression with respect to 
the coordinate Y (which is delimited by \Y\ < L y /2 due 
to the rigid confinement) so that the current becomes 

+00 

I * = l [Nf(U + V+)-N f {U + V-)] (102) 

m=0 

with V± = V(X,±L y /2). Remarkably it turns out from 
this integration that the bulk states do not formally con- 
tribute to the current. This property is a consequence of 
the continuity equation obeyed by the vortex current den- 
sity in two-dimensions. Furthermore, we sec that the op- 
posite edges (+ and -) contribute to the current with an 
opposite sign. This sign confers a chirality on the edges. 
Here, we recover from a bulk approach the importance of 
the edges for the transport properties in high magnetic 
fields that had been pointed out by several authorsi 6 ' 16 i 17 
Besides, from boundary requirements, we have to im- 
pose the condition that the current density flowing in the 
y direction vanishes at the edges 

+00 

J y = ~ n ^ + V(R)) d x V±(R) = 0. (103) 

rn — 

We obtain from this boundary condition that the edges 

Y = ±Ly/2 are at constant potentials (dxV± = 0). Let 
us take e.g. V± = ^feUn/2 where Uh is the Hall volt- 
age. Obviously, the longitudinal conductance G xx is then 
zero. Differentiating the current expression (|102[) with re- 
spect to the voltage Uh, we obtain the differential Hall 
conductance (defined by G yx = —dl/dUn) 

e 2 

G V- = ^ E 1^ + + UF + ^-)K 104 ) 

m=0 

In the linear response Uh ^ 0, we get straightforwardly 
e 2 ^ 

G yx = — 2J n F (e m - /j,), (105) 

rn— 

G xx = 0. (106) 

At zero temperature, the Hall conductance is therefore 
quantized as G yx = ne 2 /h as long as the chemical poten- 
tial e n < ji < s-n+i- This quantization is clearly associ- 
ated with a vanishing longitudinal conductance G xx = 0. 

As shown in the previous Section, the chemical poten- 
tial \x can effectively take values intermediate between 
the Landau levels e n when varying the magnetic field be- 
cause of the broadening of the density of states caused by 
the distribution of amplitudes of the impurity potential 

V in the bulk. In the non-linear regime {Uh 7^ in the 
right-hand side of the equation (|104j) ). we see that the 
Hall conductance can still be quantized at zero temper- 
ature as long as both edge states arc filled or empty, i.e. 
when the chemical potential \i lies above or below the 
pair of edges levels. 

It is worth stressing here the difference with the spec- 
tral arguments mentioned in the Introduction that are 
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FIG. 3: (color online) Hall conductance as a function of the 
inverse magnetic field for various disorder broadenings. Here 
the temperature is T = £f/200. The three disorder param- 
eters r considered are the same as in the Figs. [1] and [2] (we 
have followed the same convention for the curves). The nicest 
plateaus are obtained for the largest broadening parameter. 



usually provided to explain the absence of dissipation in 
high magnetic fields. We see from the present microscop- 
ical quantum derivation of transport properties that the 
key result leading to the vanishing of G xx and to the 
quantization of G xy is the expression of the total current 
in terms of the edge states only. Such a property has 
straightforwardly resulted from the possibility to define a 
vortex current density at high but finite magnetic fields. 
We thus deduce that the existence of electronic vortex 
states, rather than the presence of a disordered poten- 
tial, is responsible for the high-field localization proper- 
ties. In fact, the bulk disordered potential plays a role 
only at the thermodynamical level through the oscilla- 
tions of the chemical potential with the magnetic field, 
as will be developed now. 



B. Study of the width of Hall plateaus 

To study the dependence of the plateaus widths on the 
temperature or the disorder broadening, it is crucial to 
take into account the dependence of the chemical poten- 
tial on the magnetic field. For this purpose we have com- 
puted numerically the Hall conductance given by the for- 
mula (|105p including the self-consistent solution for the 
chemical potential determined from the equation (|90[) . 

For a given finite temperature, the plateaus become 
narrower when decreasing the impurity broadening pa- 
rameter, as can be clearly seen in Fig. O This paradoxal 
trend is well known experimentally and led to the con- 
clusion that the dirtiness of the sample is crucial for the 
integer quantum Hall effect. Besides, the quantum Hall 
effect does not occur in samples with strong disorder. 
The important condition to understand this limitation by 
the disorder is the condition for the high-field expansion. 



G xy [e 2 /h] 
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e F /hw c 



FIG. 4: (color online) Hall conductance as a function of 
the inverse magnetic field for various temperatures. Here 
the disorder broadening parameter F — £f/10 correspond- 
ing to a relatively flat density of states (see Fig. [T}. The 
widest Hall plateau is obtained for the lowest temperature, 
here T — £f/80. The width of the plateau shrinks with in- 
creasing temperature. The other curves correspond to the 
higher temperatures T = £f/40, T = £f/20, and T = ep/10. 

The expansion is roughly valid provided that the char- 
acteristic length of variation of the impurity potential is 
typically bigger than the magnetic length. This condi- 
tion is fulfilled only for relatively high magnetic fields 
and smooth disordered potential. We note in Fig. [3] that 
on the other hand the disordered potential has to provide 
a sufficient broadening of the density of states in order 
to get wide plateaus. 

In Fig. [5J we study the temperature dependence of the 
plateau width for a fixed disorder configuration charac- 
terized by a finite T. For high temperatures, the plateaus 
are almost invisible and we recover the classical Hall law 
for the transverse conductance. We see that the full 
plateau width is restored by decreasing the temperature. 
This trend is found for any finite T, even very small. 
Thus, a large density of states between the Landau levels 
as seemingly often the case in experiments^^ on the 
integer quantum Hall effect is not a necessary condition 
to observe wide plateaus of the Hall conductance. To ob- 
tain wide Hall plateaus, a small density of states between 
Landau levels just requires lower temperatures. Never- 
theless, at any finite temperature the nicest plateaus are 
obtained for the broadest Landau levels. 

This tendency can entirely be understood in terms of 
the location of the chemical potential when varying the 
magnetic field and thermal activation effects. For this 
purpose, it is instructive to plot the chemical potential fi 
divided by the cyclotron energy hoj c , see Fig. [5] We note 
that in the quasi-absence of states between the Landau 
levels, the ratio (j,/Huj c remains very close to the value 
n + 1/2 for a wide range in magnetic fields. The chem- 
ical potential is then separated from the energy of the 
edges contributing to the current by a very small energy 
gap, that gives rise at finite temperatures to an impor- 
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FIG. 5: (color online) The chemical potential /j, over hu> c 
as a function of the magnetic field. For a small broadening 
of the Landau levels (dotted curve), the ratio fj./hu] c remains 
close to the value n + 1/2 for a wide range in magnetic fields. 
For a large broadening (solid curve), the chemical potential 
approaches the energy (n + l/2)fru) c only for the fields cor- 
responding to ef = (n + l/2)hu c - For the correspondence 
between the different curves and the chosen parameters, see 
Figs. [HE] and H 



tant thermal activation and a disappearance of the Hall 
plateaus. On the contrary, for a significant overlap be- 
tween the Landau levels, the deviation of the Hall con- 
ductance from the quantized value by thermal activation 
effects is considerably reduced over a much wider range of 
magnetic fields, since the chemical potential approaches 
the energy of the edges more rarely. 



C. Higher order contributions to the current 

We have obtained the quantization of the Hall conduc- 
tance by considering only the contributions of the first 
order Green functions to the current formula. In fact, 
our formalism allows to derive in a systematic and rig- 
orous way higher order contributions. Our calculations 
(the details will be published elsewhere) have revealed 
that the next order terms in the current formula are pro- 
vided by the third-order Green functions. As for the first 
order calculation, it is again possible to define a vortex 
current density j(R). We find after cumbersome calcu- 
lations that the current density due to third-order pro- 
cesses in Landau-levels mixing is given by 



J 



;(3) 



+oo ]2 



h ^ 2 

rn=0 



n F (U + V) (m + l)AV + 



hlOr 



(107) 



r 



Obviously, this third-order contribution to the current 
density satisfies the continuity equation V • j*- 3 - 1 =0. By 
following the same derivation as presented for the first 
process in Landau-levels mixing, we thus deduce that the 
total current is expressed only in terms of the edge con- 
tributions also at third-order. 

In fact, from the expressions (f9"2)) . ([95]) and ([55]) . it is 
clear that it is possible to define a vortex current density 
at any order in the high-field expansion as 



j(R) = -i J2 J ^G<,,( w )(i/Ij| 



(108) 



R=R' 



Using the relation the vortex current density can 

also be written quite generally in terms of the retarded 
and advanced Green functions as 



J( R H £ / ^ in ^)(G R -G A ) v v , (z/|j| 



(109) 



R=R' 



The possibility to define a vortex current density at any 
order in the high-field expansion has very strong implica- 
tions. Provided that the continuity equation V ■ j = is 



checked at any order (a general proof of it will be the sub- 
ject of future work), we can conclude that the expression 
of the total current in terms of the edges contributions 
only is exact as long as the high-field expansion is appli- 
cable. We interpret this exact result as a fundamental 
manifestation of the robustness of the electronic vortex 
states, which gives rise to strong localization features in 
the quantum transport properties. 



VI. DISCUSSION AND CONCLUSION 

In this Section, we would like to discuss some impor- 
tant issues concerning the quantum Hall effect in the light 
of the theory developed in the former Sections. 

A. Longitudinal transport 

When the Hall resistance R xy jumps from one Hall 
plateau to the adjacent one, peaks in the longitudinal 
resistance R xx appeari^ A standard picture to explain 
such a behavior of the longitudinal conductance is that 
at zero temperature the edge of the electron gas generally 
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will follow some trajectory near the physical edge which is 
only on average parallel to it when the Hall conductance 
is in the plateau region. As the step in the Hall conduc- 
tance is approached, the equipotcntials have excursions 
ever deeper into the bulk of the sample, and at a critical 
value of the Fermi energy, there is an equipotential con- 
tour that connects the two edges. This is supposed to be 
the origin of the peak in the diagonal conductance. 

Within our present one-body approach, we have found 
that the longitudinal conductance is zero. The prob- 
lem with this zero longitudinal conductance has resulted 
straightforwardly from the Eq. Q103P that expresses the 
absence of current flow accross the outer surfaces. To 
our opinion, this physical boundary condition can not 
be questioned. Indeed, it is the essential physical condi- 
tion by which the presence of boundaries is accounted in 
the purely classical derivation of the Hall effect (see e.g. 
Ref . lill ) . In a formalism making use of wavefunctions the 
presence of boundaries will be introduced with the condi- 
tion that the wavefunction has to vanish at the edges. In 
a Green functions formalism the edges are usually taken 
into account via boundary conditions on the current (see 
e.g. Ref. concerning the theory of mesoscopic super- 
conductivity with quasi-classical Green functions). This 
is precisely the approach chosen in our theory, with the 
peculiarity here that the boundary conditions are applied 
on the vortex current instead of the electronic current. 
Such an approach assumes that there exists some degree 
of freedom characterizing the Green functions that is only 
completely determined through the boundary conditions. 
Coming back to the Dyson equation (|6^|) this indetermi- 
nacy of the Green functions is however hidden in the 
electrostatic part of the potential V. 

We note that the transition between Hall plateaus 
falls within the description of quantum phase transi- 
tions. Clearly, the transition between the Hall plateaus 
is highly sensitive to any fluctuations present in the sys- 
tem. We remark that the explanation of the longitu- 
dinal transport by fluctuations of the electronic density 
has also been proposed in several worksi 43 i 44 ' 45 i 46 Recent 
experimental 4 ! and theoretical 4 ^ works have focused on 
the phenomenological relationship between the diagonal 
and the transverse resistances observed experimentally, 
namely R xx oc B x dR xy /dB. From these works it turns 
out 4 ! that all R xx features could be understood in terms 
of electron density inhomogencitics (at least in the regime 
of the fractional quantum Hall effect). Our theory has 
formulated the problem of transport in terms of vortex 
motion. However, inhomogeneities of the vortex density 
have not been considered in our calculations. 

The above discussion clearly calls for a proper treat- 
ment of the electrostatic problem. By including electron- 
electron interactions in our formalism, the inhomo- 
geneities of the vortex densities translate into a self- 
consistent determination of the local electrostatic poten- 
tial (e.g. in a Hartrcc-Fock approximation) for which the 
boundary conditions on the current have to be imposed. 
This can be understood as spatial variations of the local 



chemical potential /i(R) at local equilibrium. Clearly, 
according to expression (|105[) . such fluctuations of the 
chemical potential at the edges create, when approach- 
ing the transition points between two Hall plateaus, parts 
of the samples characterized by a different quantization 
of the Hall conductance. This automatically implies that 
the edges are no more equipotential lines, so that there 
exists a longitudinal voltage difference. A detailed and 
quantitative study of this problem within our formalism 
is beyond the scope of the present paper and is therefore 
postponed for future work. 



B. Breakdown of the integer quantum Hall effect 

The quantum Hall effect is known^i 3 - to break down 
at low magnetic fields and for high current intensities. 
These limitations can be captured qualitatively within 
the condition of applicability of the high magnetic fields 
expansion theory. Indeed, the quantum Hall effect regime 
holds roughly as long as the potential V (which includes 
the electrostatic potential and the impurity potential) is 
smooth at the scale of the magnetic length, i.e. Ib <C A 
where A is the typical length scale for the potential. An 
increase of the current intensity induces unavoidably an 
increase of the net gradient of the electrostatic poten- 
tial, and thus a reduction of the characteristic length 
scale A of the potential variations. Equivalently, the in- 
equality Ib <C A breaks down in low magnetic fields since 
Ib oc B~ x / 2 diverges when B — ► 0. The quantitative de- 
termination of the critical values of current and of mag- 
netic field seems arduous since it supposes the accurate 
determination of the radius of convergence of the func- 
tional series an d of the high-field expansion of the 
Green functions. 



C. Relationship between the quantum Hall effect 
and the Shubnikov-de Haas effect 

It is worth emphasizing that the low magnetic field 
regime is usually characterized by another quantum ef- 
fect, namely the Shubnikov-de Haas oscillations of the 
longitudinal resistance R xx , while the high magnetic field 
regime corresponds to the quantum Hall effect. In the 
Shubnikov-de Haas regime, R xx presents minima between 
the Landau levels that become deeper and deeper when 
increasing the strength of the magnetic field, whereas 
there is a preformation of rough steps in the R xy fea- 
tures. In the quantum Hall effect regime, R xx drops to 
zero. The passage from the Shubnikov-de Haas regime to 
the quantum Hall regime is relatively smooth. It suggests 
the existence of a physical relationship between these two 
quantum effects. 

The theory of the Shubnikov-de Haas effect for the two- 
dimensional electron gas in the presence of a disordered 
potential has been developed a long time ago within the 
Green functions formalism and the standard impurity av- 
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craging procedure ! i 50 The oscillations of R xx are usu- 
ally understood in terms of the oscillations of the density 
of states. Within this picture, the minima of R xx are re- 
lated to the minima of the density of states. Theoretical 
derivation n 37 ' ' 52 ' 53 ' 54 of these oscillations making use of 
semi-classical Green functions provide a deeper interpre- 
tation concerning the physical origin of the Shubnikov- 
dc Haas effect. Within this semi-classical approach, the 
oscillations of R xx and of R xy are interpreted as interfer- 
ence effects due to a fraction of electrons following clas- 
sical periodic orbits, in close relation to the Gutzwiller's 
trace formula^ for the density of states^ This fraction 
P = cxp(—2ttR c /1) (where I is the mean free path) of 
electrons participating to magnetic quantum effects in 
low fields^ 7 - corresponds to the classical probability that 
the electron will complete a cyclotron orbit without hav- 
ing suffered a collision.— Note that the disorder tends to 
reduce this probability and thus destroys these magnetic 
interference effects. This is in complete contrast with the 
weak localization effects that are induced by the disorder 
and are suppressed by the application of a weak mag- 
netic field (see e.g. Ref. Since the cyclotron radius 
diminishes with the magnetic field, the fraction of circling 
electrons increases and the minima of R xx become more 
and more pronounced. In the limit of infinite magnetic 
field, we expect from the semi-classical picture that all 
the electrons participate to these magnetic interference 
effects. 

In the present work, we have shown that the quanti- 
zation of the orbital degree of freedom confers a robust- 
ness on the periodic cyclotron motion, so that the elec- 
tronic states have a probability higher than the classical 
one (given above) of returning to a same given position. 
The ineffective energy exchange between the orbital de- 
gree of freedom and the vortex degree of freedom occurs 
consequently in any local region of smooth potential in 
high but finite magnetic fields, i.e. not only in the semi- 
classical limit B = 00. This property resulting from the 
topological character of the circulation quantum number 
m is illustrated by the possibility to define a vortex cur- 
rent density j (R) at any order of the high magnetic fields 
gradient expansion. The electronic system enters into the 
quantum Hall regime when all the electronic vortex states 
contribute to the magnetic interference effects (strong lo- 
calization). Accordingly, the Shubnikov-de Haas effect 
can effectively be seen as a precursor to the quantum 
Hall effect. 



D. Conclusion 

The insensitivity to the type of impurities and the re- 
productivity of the experimental results independently of 
the type of host material has suggested that the quantum 
Hall effect is due to a fundamental principle. In the first 



half of the paper, we have shown that the quantization of 
the energy spectrum into Landau levels is associated with 
the topological condition of the singlc-valucdncss of the 
single-particle wavefunction. Electronic states can thus 
be characterized by vortex states whose nature comes 
out in the presence of a perturbation (for instance a dis- 
ordered potential). As developed in the second half of 
the paper, the robustness of the vortex states is respon- 
sible for the remarkable quantum transport properties of 
a two-dimensional electron gas in high finite magnetic 
fields, namely the vanishing of the longitudinal conduc- 
tance and the quantization of the Hall conductance. 

By reformulating the high-field expansion theory us- 
ing a vortex states representation, we have rigourously 
derived within the formalism of real-time Green func- 
tions a quantum mechanical expression for the current. 
Within our vortex representation, a current is obviously 
associated with a (hydrodynamic) vortex flow. We have 
shown that this quantum mechanical property is an ex- 
act result that holds as long as the high-field expansion is 
applicable. We have demonstrated explicitly from a bulk 
derivation up to third order processes in Landau-levels 
mixing that the total current can be expressed entirely 
in terms of edge contributions. This remarkable prop- 
erty together with boundary requirements leads directly 
to the absence of dissipation and to the quantization of 
the Hall conductance. 

We have emphasized that the density of states is not 
necessarily small between the Landau levels in the vortex 
regime. By computing the oscillations of the chemical 
potential with the magnetic field, we have pointed out 
that at any finite temperature the widest widths for the 
Hall plateaus are precisely obtained in the case of the 
broadest Landau levels, in agreement with experimental 
observations. 

The quantum mechanical picture and the formalism 
developed in the paper are appealing to explain several 
issues concerning the quantum Hall effect. We have al- 
ready pointed out qualitatively the relationship between 
the Shubnikov-de Haas regime in low magnetic fields and 
the quantum Hall regime in higher magnetic fields. More- 
over, as also underlined, our formalism allows a system- 
atic calculation of higher order contributions to the cur- 
rent. Finally, we hope that it will be possible to gener- 
alize the presented approach of the quantum transport 
to incorporate the interaction effects between the elec- 
trons that appear to be crucial to capture the dissipative 
longitudinal transport. 
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APPENDIX A: OVERLAP 

In this Appendix, we calculate the overlap between the two vortex states |m, R) and |m',R') 

(m,R|m',R) = J d 2 r% lR (r)* m ,, R ,(r). (Al) 
Making first the change of variables u = r — (R + R')/2 and using the cartesian coordinates u = (x,y), we obtain 



. Rim ,R) — C m C m ' 



d 2 u 

I 2 

' B 



(x + d x ) -i(y + d y ) 



V2h 



(x - d x ) + i(y - d v ) 



-1 m (x-idyy + ( y +id x y 



V2h 



(R-R^) 2 -2iz-(R.xR / ) 



XC 



(A2) 



with d = (R' — R)/2. We then make the shift of the complex arguments id y and id x in the x and y integrations 
respectively. Because the functions are analytic it is possible to deform the contour to the real axis so that 



(m,RK,R') = CX 
The remaining integrations yield the result 

(m, R|m',R') = cxp 
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' x - iy 
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(A3) 



(A4) 



The overlap contains a phase factor resulting from the interference effects between the states \m, R) and \m, R'). 
The basis of vortex states becomes orthogonal with respect to the quantum number R only in the limit of infinite 
magnetic fields l B — ► 0. 

In fact, from a superposition of the states |m, R) with different vortex positions, it is possible to construct orthogonal 
sets of eigenstates at finite magnetic fields. For example, the states 



dY^ r 



R 



(A5) 



obtained from the superposition of states \m, R) (expressed in the symmetrical gauge) with their vortex positions on 
the axis Y form an orthogonal set. Inserting the expression (|34[) of the vortex states we perform the integration over 
Y in Eq. (|A5[) using the identity 



/+oo 
X m cxp [—(A - (3) 2 ] 
-OO 



{2i) n 



H m {if3) 



(A6) 



where H m is the Hcrmite polynomial of degree m. After calculation, we find that the states (|A5|) can be expressed as 

(x-X/2) 2 ] u fx-X/2- 



*m,x(r) 



KV2) 



m—l 



_ iyX 

e 21 ~b exp 



4Z 2 



h 



(A7) 



which is nothing else than the expression of the (non normalized) Landau states (J5J written in the symmetrical gauge. 
From the comparison between the expressions (|A7[) and ([5]), we can identify p = —X/2l 2 B and xq = A/2. 



APPENDIX B: DERIVATION OF THE KELDYSH GREEN FUNCTION 



In this Appendix we derive the equations obeyed by the Green functions in the Keldysh space within the vortex 
states representation. The quantum kinetic equation for the Keldysh Green function is then solved within the high 
magnetic field expansion theory. 
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1. Dyson equation in the Keldysh space 

The Green functions in presence of the perturbation V are obtained from the Dyson equation 

G = G + G VG (Bl) 
where G designates the 2x2 matrix in the Keldysh space defined as 

G R G K 



G ={0 G A )■ ^ 

Here V may contain a disordered potential and an electrostatic potential. Inserting the quantum numbers \v) = \m, R) 
and considering that V is local in time, we have 

G(i>i,ti;i>2,t2) =G (i>i,ti;ifii,t2)+ ^ fdt 3 Go(vi,ti;v 3 ,t 3 )V V3>Ui G(v4,t 3 -,V2,h)- (B3) 



^3,^4 



By applying the time derivative to the Dyson equation, we thus obtain the system of equations for the three Green 
functions 

i-^--Zm 1 )G R ' A (v 1 ,t 1 ;v 2 ,t 2 ) = (^i|^)5(ii-t 2 )+^V; i ,„3G i? ^( l /3,ii;^2,i2), (B4) 

h'J- - £ mi J G K (v 1 ,t 1 ;v2,t 2 ) = ^2v Vl ,„ 3 G K {v 3 ,t v ,V2,t2), (B5) 

where we have used the equation of motion for Go 

d \ ~ 

i dT~^ mi ) ^o^ 1 '* 1 ^ 2 '* 2 ) = ( v iW2)5(ti-t 2 )I (B6) 

with / the unit matrix in the Keldysh space. From the other Dyson equation (G = Go + GV^Go) and by applying the 
derivation with respect to the time t 2 , we get another equation for the Keldysh Green function 



(B7) 



It is direct to see that the retarded and advanced Green functions depend only on the time difference r = t\ — t 2 . 
Fourier transforming with respect to t\ — t 2l we get the equations 



(w - U ± i8) G^» = ( Vl \v2) + V»u»*Gu£M- ( B8 ) 



As for the Keldysh component, we take the sum and the difference of Eqs. (|B5p ~ (|B7|) . introduce the relative time 
t = t\ — t 2 and the time T = (t\ + t 2 )/2, and Fourier transform with respect to r, which yields 

(2u-i$ mi +U 2 }) G^(r,w) = ^[V^ lVa G^(T,w) + G^(T,u)V^ lVa ] , (B9) 
(^-[U-fm,]) G« uV2 {T,w) = Y, [Vu l ,» 3 G«, V2 (T,u;)-G« uV3 (T,u;)V V3 , l/2 \ . (BIO) 

2. Derivation of Keldysh Green function in the high magnetic field expansion 

As for the retarded and advanced Green functions in Sec. III.C, the Keldysh Green function is searched under the 
form of an expansion in terms of the order of (7b /A) 9 . The leading order contribution obeys the equation 
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+UD Gg£(T, W ) = X; , (Bll) 



i^=- [U -U]) G^£(T,u;) . (B12) 



<9T 

For the second terms in the right-hand side of these two latter equations, we can take the potential V at the point 
R2 instead of R3 since at leading order we disregard correcting terms of the order of Zs/A. Therefore we get 

(2cj - [£ mi + U 2 ] - V(Ri) - V(R 2 )) G<$K{T,u) = 0, (B13) 
(i-^ - [£ mi - £ m J - ^(Ri) + ^(R 2 )) GV>£{T,w) = 0. (B14) 

Disregarding higher order correction terms, this system for the zero order Kcldysh Green function can in fact be 
written as 

(2u - [£ mi + C™ 2 ] - 2V(Ri)) G%>£ (T, w) = 0, (B15) 

1-^ - [u ~ uS) gMJZP, w) = 0. (Bi6) 

The equation (|B15P for the zero order Keldysh Green function yields straightforwardly the following spectral depen- 
dence 

G<°> * (T, w) « S (cj - [U, + - V(Ri)) . (B17) 

The differential equation (|B16j) is readily integrated: 

Gg>£ (T, cj) = G„ 1)1/2 (w) exp (i[U2 - UiYO ( B1 8) 

where C P1M2 (u>) has to be determined from an initial condition. Taking into account that the perturbation V (including 
the external field E) is switched on adiabatically with the equilibrium initial condition at T = — 00 

G£)£(-oo,w) = Gff(w) = -2i7rtanh(a;/2T) (y x \v 2 ) 8{lo - f mi ), (B19) 

we obtain 

Gg>£(7>) = -2i7rtanh(o;/2T) ( Vl \u 2 ) 5 (w - f mi - f (R x )) = tanh(w/2T) (g<°>* - G^^) . (B20) 

Therefore, at leading order the Kcldysh Green function is still diagonal with respect to the quantum numbers m. 

We shall determine the higher order Keldysh Green function G^ K within the stationary regime. Then, the first 
order Keldysh Green function is determined for example from the equation resulting from the sum of Eq. (|B10[) and 
(|B9| which yields up to first order in Ib / A 

(u-U - V(Rx)) G£>£(cj) = ^KW 3 (R 1 |R 3 )G(0)^(cj). (B21) 

Inserting the zero order Green function and performing the sum over 1/3, we get 

G (1 > K = tanh(w/2T) (g (1 > r - G (1) A ) . (B22) 

We see that the structure of the Dyson equation for the Keldysh Green function of orders higher than one in the 
stationary regime is similar to the structure of the Dyson equation for the retarded and advanced Green functions 
at the same order. Consequently, we can deduce that the proportionality (|B22[) between the Keldysh Green function 
and the retarded and advanced Green functions holds at any order of the high-magnetic field expansion theory. 
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APPENDIX C: MATRIX ELEMENTS OF THE CURRENT DENSITY OPERATOR 



In this Appendix, we derive the matrix elements of the current density operator in the vortex states representation. 
In the space representation (|r)), the current density operator is written as 



Adopting the symmetrical gauge 



j(r) = 



A = l(-y 

2 



he l 
H27| 



we express the matrix elements of the current density operator as 
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i\frnJ 6 m ,m'-i - i\fm.8 m , m i + i 
r m' S m ,m'-i + \/rri8 m ,m'+i 



(CI) 

(C2) 

(C3) 
(C4) 
(C5) 



Note that these matrix elements of the current density operator in the vortex states basis couple adjacent Landau 
levels exactly as in the basis of Landau states. 
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